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Abstract. Fhat the RR Lyrae star period - K-band luminosity relation is a promising tool as a distance indicator in the Milky 
Way and Local Group of Galaxies is apparent on observational and theoretical grounds in the literature. Less clear is the sensi- 
tivity of the relation, and consequently the physics of horizontal branch stars, to differences in stellar environment. In this paper, 
the first measurement of the (fundamental) period - K-band luminosity relation for the central region of a globular cluster is 
presented. It is based on a sample of seven RR Lyrae stars imaged with adaptive optics. In addition, the relation for the outer re- 
gion has been reanalyzed, and is found to be in good agreement with both the previous estimate by Longmore et al. (1990), and 
with the inner region relation, especially when irregular and double-mode RR Lyrae stars are excluded. Importantly, there is no 
difference between the slope of the inner and outer region relation within measurement uncertainties, suggesting no difference 
in evolutionary state (luminosity). Faking the M3 distance modulus as 15.0±0.07 mag, the period-absolute K-band magnitude 
relation derived by linear least squares fitting is: 

Mk = —0.96 (±0.10) — 2.42 (±0.16) LogP for the inner region. Excluding irregular variable stars, the outer region rela- 
tion is: Mk = —1.07 (±0.10) — 2.38 (±0.15) LogP . Fhis good agreement provides further strong support for the near-IR 
period-luminosity relation as a distance indicator. 

Key words. Stars: variables: RR Lyr - Stars: distances - Stars: horizontal branch - Stars: imaging - Instrumentation: adaptive 
optics 



1. Introduction 

For more than a century, RR Lyrae stars, in addition to classical 
Cepheids, have played an important role as primary distance 
indicators in the Milky Way and the Local Group of galax- 
ies. Fhe importance of RR Lyrae stars stems primarily from 
their brightness, being presently observable efficiently out to 
distances of several Mpc; they have characteristic light curves, 
easily measured periods, a narrow magnitude range, and a high 
frequency in many globular clusters, but their intrinsic luminos- 
ity is uncertain, both empirically and theoretically (see Caputo 
et al. 2000, and reference therein). Indeed, the dependence of 
intrinsic luminosity on metal content is uncertain: metal-rich 
and -poor RR Lyrae stars appear to follow different linear My 
- [Fe/H] relationships (Caputo et al. 2000). 

In contrast, near-infrared observations of RR Lyrae stars 
hold several distinct advantages over optical investigations: 
they have a significantly weaker dependence on [Fe/H] (Bono 
et al. 2001); the K-band amplitudes of RR Lyrae variables are 
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smaller than their V-band values; and their near-IR light curves 
are more symmetrical. All these points lead toward more easily 
measured and accurately determined absolute magnitudes, and 
importantly, with a much reduced sensitivity to reddening. The 
relationship between the periods and K-band magnitudes of RR 
Lyrae stars has been measured in the pioneering semi-empirical 
work of Longmore et al. (1990) who found a tight relationship 
for several galactic globular clusters. Bono et al. (2001) pro- 
vided the theoretical support for these findings. They calculated 
the dependence of the relation on metallicity, and determined 
the intrinsic dispersion in the three relevant parameters, namely 
period, luminosity and chemical composition. Confirming the 
reliability of the relation further, Bono et al. (2001) applied the 
PLk relation together with predicted mass and luminosity from 
current evolutionary models, to forecast the parallax of RR Lyr 
and found a good agreement with the HST result (Benedict 
2002). 

As quoted in the literature before, horizontal branch (HB) 
stars could be regarded has having crucial importance in our 
understanding of several unresolved astrophysical problems 
such as the second-parameter problem, the UV-upturn in ellip- 
tical galaxies (Ferguson 1999), as well as the dependence of the 
absolute V-band magnitude My of RR Lyrae stars on chemical 
composition. 
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A relevant key issue is abundance variations in stars. 
Abundance variations (typically of C, N, and O) have been 
found to occur along the whole stellar evolutionary sequence 
from below the turn-off of the main sequence to the red giant 
branch (RGB) in several globular clusters. It is hard to list all 
the significant reviews in this field, but Kraft (1994), Cannon 
et al. (1998), Briley et al. (2001) and references therein pro- 
vide an overview of current thoughts on the subject for giant 
stars, while on the topic of main sequence stars, useful refer- 
ences can be found in Harbeck et al. (2003). RR Lyrae star 
luminosities, colours, and periods could be affected by inho- 
mogeneities in the initial chemical composition of the primor- 
dial gas, internal stellar mixing processes, mass loss, and pos- 
sibly multiplicity in globular clusters, and the cumulative ef- 
fect of such variations may reveal themselves in precise multi- 
wavelength period-luminosity relation studies. It is fair to say 
that a database of precise, consistent measurements for a large 
number of RR Lyrae stars covering a range of environments 
(in globular clusters, the field, different galaxies, etc.) will be 
important for useful future tests of the intrinsic dispersion pre- 
dicted by Bono et al. (2001). 

In this paper, the period-K-band luminosity relation for the 
inner region (r < 20") of the globular cluster M3 is presented. 
As a close agreement with the outer region relation (r > 50") 
is generally expected to exist but has not been shown, the in- 
ner and outer region relations are compared. The motivation 
has been to provide further support for the near-IR period- 
luminosity relation as a distance indicator. This work has been 
made possible by K-band adaptive optics - assisted observa- 
tions. 

2. Observations and data analysis 

The data consists of a series of K-band frames taken on April 
19 th , 21 st 2000 using natural guide star adaptive optics (AO) 
imaging with ALFA (Kasper et al. 2000) on the 3.5 m telescope 
at the German-Spanish Astronomical Centre on Calar Alto. 
The science camera was Omega Cass with a plate scale of 0.08" 
pixel -1 giving a field of view (FoV) of 80" x 80". Dithering 
of source frames during exposures was not performed, but 
there were occasional changes of guide star between exposures. 
Observing overheads were of the order of a few xl0%, and 
there were target-of-opportunity override observations at the 
beginning of each night. Conditions were photometric on the 
19 th and some light cloud was present on the 21 st . The K-band 
source details are given in Tabled On the 21 st some J-, and H- 
band frames were taken, but are omitted as they are not within 
the scope of this study. The Shack-hartmann wavefront sensor 
(WFS) FoV was set to 3" to block out neighbouring stars, and a 
signal threshold level was used to ignore faint stars in the WFS 
FoV. There were two suitable AO guide (or reference) stars, 
namely SI (my ~ 12.8, mx = 9.3) and S2 (my ~ 13.2, mx = 
10.2), marked in Fig. ^ A consequence of the AO correction 
is that the best correction occurs in the direction of the guide 
star, becoming worse at off-axis angles. There is no obvious 
evidence of systematic residuals due to radial PSF changes in 
the star-subtracted frames, though the central PSF flux is more 
pronounced toward the guide star. Through comparison with 



K-band data from Davidge & Grundahl, in preparation, which 
is reported in the next section, the effect on photometric mea- 
surements is up to a few percent at most, similar to reports 
by Davidge (2000) and references therein. The point spread 
function (PSF) FWHM varied between frames from about 0.3 
to 0.6" (see Table Q due to varying seeing, and signal-to- 
noise. The globular cluster M3 (NGC 5272; d ga i ac t ccntric = 
12.2 kpc; Harris 1996) was selected because of its plentiful sup- 
ply of RR Lyrae stars (Corwin et al., in preparation; Corwin 
et al. 2001; Clement 1997; Sawyer-Hogg 1973). Several po- 
tential guide stars are within the central few arc-minutes, and 
there is a relatively low number density of stars compared to 
M15, which presently is the only other cluster with both a well 
studied, and plentiful central RR Lyrae star population (Butler 
et al. 1998; referred to hereafter as B98). Additionally, pub- 
lic Hubble Space Telescope (HST) integrations taken in 1995 
April 25 are available and were retrieved from the STScI sci- 
ence archive. Calibration details are given by Bagget et al. 
(2002) and Holtzman et al. (1995). The observations consist of 
two 70s and two 3s Wide Field Planetary Camera 2 (WFPC2) 
integrations with the F336W and F555W filters respectively 
(referred to here as U336 and V555). They are listed in Tabled 
The U336 data was used to provide accurate stellar astrometry, 
and both filters were required for high signal-to-noise giant star 
magnitudes. 

In the U336 frame, the core of M3 is roughly centred on 
the Planetary Camera (PCI) chip which has the highest angu- 
lar resolution (0.045" pixel -1 ). It was thus decided to concen- 
trate on the PCI frame only. The exact locations of stars have 
been obtained by reducing the Planetary Camera (PCI) frame 
following the procedure developed by B98. The resulting stan- 
dard deviation (rms) of the astrometry residuals from a high or- 
der U- to V-band coordinate transformation are of the order of 
0.0025" in both axes. By using this astrometry the positions of 
stars in the K-band field have been accurately located. For the 
present application, the simple charge transfer efficiency ramp 
correction of Holtzmann et al. ( 1995 1 is adequate and has been 
applied. The U336 and V555 have been zero point calibrated ac- 
cording to the Holtzman et al. i 19951 1 recipe (their Eq. (8) and 
Table 7). 

For the near-infrared data reduction standard IRAF 1 rou- 
tines and scripts were used. The basic image frames (typi- 
cally 60s) were taken with the cluster centre positioned on 
one of the CCD quadrants, and a standard near-infrared reduc- 
tion scheme of sky-subtraction and flat-fielding was followed. 
However, the measured sky frame illumination is non-uniform, 
varying across each camera quadrant, predominately along the 
readout direction. In contrast, the dome frames have a flatter 
appearance. The effect is such that the drop-off in illumina- 
tion intensity is up to 20% in each quadrant. The cause is not 
clear. For correction of this effect, the parameters involved are 
the number of counts per pixel from the source S, the thermal 
signal T, and the observed signal O which is related to the nor- 



1 IRAF is distributed by the National Astronomy Observatory, 
which is operated by the Associations of Universities for Research 
in Astronomy, Inc., under the cooperative agreement with National 
Science Foundation. 
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a Measured at the centre of the PCI region, about 20" from the guide star(s); Uncertainty in spatial resolution is ~ ±0.05". 




Fig. 1. Composite/cropped frame showing the central (K-band) 
field of M3 with the two AO guide stars, SI and S2., and rele- 
vant RR Lyrae stars marked. North (N) and East (E) directions 
are included. 



malized flat-field response (or relative quantum efficiency) of 
the imaging system, k, by 



0(x,y) = k(x,y)xS + T(x,y) 



(1) 



The standard iterative procedure of fitting a low order surface 
function to the median filtered quadrants was applied to sky and 
dome flat-fields to determine the conversion image k s k y / k<j ome 
(where Sdomc.sky = 1-0). For each source image, the final flat- 
field was created by taking the inverse-normalized lights 'on' 
and lights 'off dome frames in turn, then multiplying each by 
the appropriate conversion image, and finally taking the dif- 
ference between the two of them. Importantly, this process re- 
moves the average bias and dark count in each pixel. 



The next issue is sky subtraction. As K-band sky emission 
varies over time, sky frames taken close in time (less than about 
two minutes) to each object frame, which might otherwise 
swamp the light from the faintest stars, need to be subtracted. 
Sky frames were built by stacking image frames taken at dif- 
ferent pointings at 300" from the cluster centre and dithered 
by 10-15". Incomplete star removal in the bottom-right chip 
of some median sky frames occurred due to the presence of a 
few close stars; but these stars are outside the PCI region and 
thus have no effect on the results presented in this paper. The 
best fitting function [kxS] s k y (S s k y = 1-0) was determined for 
each stacked sky frame; subtraction of it from the sky frame 
itself yields the thermal component. As it is important to re- 
move most of the night sky flux, the sky component (the best 
surface fit to the appropriate sky frames) was scaled slightly for 
most science frames, typically by about 1%. The thermal com- 
ponent, which is about 1-2% of the sky flux was then added. 
For hot/bad pixel removal, the flat-field was multiplied by a hot 
pixel mask derived from the flat-fielded sky frame. Next, for 
sky background removal, the flux-weighted centroid of each 
source frame was determined and the frames were shifted to 
a common centroid. This stacked source frame and the stacked 
sky frame, both of equal integration time, were then flat-fielded 
and the difference was taken. 

Following the calibration of all the K-band frames, eleven 
frames with sub- half arc second spatial resolution were se- 
lected for signal-to-noise/crowding reasons, the flux-weighted 
centroid of each frame was determined, and all were shifted to a 
common centroid. This was done three times using a different 
set of frames at each turn to allow three high signal-to-noise 
estimates of the K-band magnitudes. In this way, a bootstrap 
estimate of magnitude uncertainty can be obtained later after 
determining the photometry. To test the internal photometric 
precision later using the brightest stars the individual frames 
were also reduced. 

It is possible that analysis of the non-PC 1 portion of the full 
K-band field in Fig.^may reveal some evidence for red giant 
branch variable stars that could be considered in the context of 
giant star evolution; but such a study is outside the scope of this 
paper. 
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2. 1. Photometry 

The photometric reductions have been carried out using the 
DAOPHOT (Stetson 1987) photometry package. The proce- 
dure used to derive U336 and V555 instrumental magnitudes, 
is taken from that described by B98. The V555 astrometry is 
estimated to be complete down to V555 = 17 (K < 15 for (V- 
K) g iants < 2) based on insertion of artificial stars into the V 555 
frame and repeating the data reduction. 

PSF-fitting was performed using the routine ALLSTAR 
(Stetson & Harris 1988), part of the DAOPHOT package. Each 
K-band frame has been considered as two sections overlap- 
ping by 20", and the following procedure for each section, 
taken partly from B98, has been carried out: The initial PSF 
was obtained using standard recursive techniques; the star list 
at this point consisted only of stars detected by DAOFIND, 
which typically found about 250-300 in the PCI region. This 
allowed a high-order coordinate transformation between the K- 
band frames and the HST/PC1 field, thereby accounting for rel- 
ative geometrical distortion near known RR Lyrae stars in the 
field. 60-90 stars were selected from each section. After apply- 
ing this transformation, one for each section, to the HST/PC1 
star list, the rms positional accuracy was of the order of 0.02" 
in both axes. Then, a background-subtracted image and its PSF 
were fed into the PLUCY (Hook et al. 1994 1 restoration task in 
IRAF. Although not photometrically reliable, the data from the 
restoration task provided good initial magnitude estimates for 
the brightest stars in the field. Together with these magnitudes, 
accurate coordinates, sky background estimates and a unique 
ID number for each star, a final PSF fitting could be performed. 
Next, U336 2 background estimates were scaled up appropri- 
ately, followed by re-making the PSF, and transformation of the 
list fitted by ALLSTAR using the new star IDs, in groups of 500 
stars; the sky background is re-calculated for each star individ- 
ually rather taking a group estimate as the background varies 
significantly across the field. Finally, the photometry from each 
section was tied to the same internal zero point. This step was 
performed by taking a calibration of a sequence of stars in the 
range K < 17 against their counterparts in an average dataset. 
Firstly, the histogram of individual offsets, typically forming a 
normal-like distribution, was created. Then, in an iterative pro- 
cess, histogram outliers were ignored until two different esti- 
mates of the zero point offset, the weighted-histogram average 
and median value, differed by less than 0.01 mag. In this way, 
the calibration of the K-band magnitude zero point is on a ro- 
bust statistical footing, insensitive to large amplitude variable 
stars and photometry outliers. 

In the context of photometric precision, it is worth mention- 
ing that a space-invariant empirical PSF made using 2 - 4 of the 
brightest stars, and modeled on a moffat elliptical function of 
index 1 .5 works best because the background light due to the 
faint unresolved stellar population is distributed in an inhomo- 
geneous way. As there is a pointing difference of about 12" for 
some frames, it was possible to test the photometric precision 

2 The light from the undetected stellar population is dominated by 
light from upper main sequence stars in a narrow colour range; the 
deviation from the mean colour for the bulk of these stars over the 
field is expected to be up to a few tenths of a magnitude. 



along both axes for each individual frame; for this, the way 
in which magnitudes differ from their time averaged/intensity 
weighted average was examined. No obvious indication of a 
trend with field position was found. 

For validation of the quality of the instrumental magnitudes 
a sequence of stars in the range K < 17 was calibrated against 
their counterparts in Canadian-France-Hawaii Telescope AO 
data from Davidge & Grundahl (2002, in preparation; referred 
to here as DG). Davidge & Courteau (1999) reported an ab- 
solute zero point error of ±0.009 mag for their K-band mag- 
nitudes which came from CFHT AO observations; but conser- 
vatively, ok = 0.05 mag is assumed for the DG data. The dif- 
ference between the intensity weighted average Kpci and the 
DG data are compared in Fig. |2j and known RR Lyrae stars 
from Corwin et al. (2002, in preparation) are marked with open 
circles. Apart from the outlying RR Lyrae star which may be 
a false measurement due to stellar crowding, there is a good 
agreement. In Fig.[3]the rms scatter is plotted as a function of 
K for stars in the PCI region with at least two (out of three) 
data points (bottom); and (top) for stars in common with the 
Kdg data for which at least two Kpci data points have been 
measured. It can be seen that the rms value in the RR Lyrae 
star magnitude range is about 0.05 mag, caused by both intrin- 
sic variability and photometric error. It is useful to note that in- 
trinsic K-band amplitudes of RR Lyrae stars of type RRab are 
up to a few tenths of a magnitude, and smaller for RRc-type 
variables. 

For zero point validation, the Cohen et al. (1978; referred 
to hereafter as C78) giant branch data was considered. The 
V-K colour of their M3 giant branch ridgeline at Mk = -3.0 
(<tv-k=±0.03 mag) was checked and compared with the value 
found on the Ferraro et al. (2002; referred to hereafter as F00) 
ridgeline for globular clusters with a similar metallicity to M3; 
[Fe/H] = - 1.34 dex from Carretta & Gratton (1997) is adopted. 
A good agreement was found to within ±0.05 mag. Next, in or- 
der to convert the F00 ridgeline Mk to K-band magnitudes, a 
distance modulus of 15.00±0.07mag (Caputo et al. 2000) (see 
Sect. 13. II later for a discussion) was assumed and the V/K data 
was overlaid with the C78 V/K giant star data on the (V, V-K) 
colour magnitude diagram in Fig.|4] A solid line fit to the C78 
data (excluding the outliers (asterisks)) is included. The slight 
mismatch in ridgelines close to the tip of the RGB arises be- 
cause the corresponding metallicity of the F00 line is only sim- 
ilar but not equal to that of M3. The lower end of the RGB is 
much less metal-sensitive, and so, as confirmed above, it is safe 
to compare the ridgelines at two or more magnitudes fainter 
than the RGB Tip. Taking the C78 data as the reference dataset, 
the offset between the C78 and F00 lines suggests that the lower 
end of the C78 ridgeline comprises asymptotic giant branch 
(AGB) stars at V > 13.5. Next, the (V 555 , V 5 5 5 -Kpci) colour 
magnitude diagram for stars with <tk < 0.03 mag and cry. 555 < 
0.05 mag (the V-band zero point error is ±0.05 mag) was over- 
laid with the Kpci magnitudes shifted by -0.1±0.05mag to 
match the lower end of the RGB fiducial ridgeline for [Fe/H]=- 
1.34 dex, (dot-dashed/thick) which was deduced from the F00 
(V-K)-[Fe/H] relations. 
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Fig. 2. Difference between the intensity weighted ensemble av- 
erage K-band photometry from this work and DG as a function 
of K-band magnitude. The RR Lyrae stars in common between 
the two data sets are marked with open circles. The circled out- 
lier may be a false measurement due to stellar crowding. 
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Fig. 3. ck is plotted as a function of mean Kpci magnitude. 
Top: Stars in common with the DG K-band data with at least 
two data points. The five known RR Lyrae stars in common 
are marked with open circles in both plots at K~15; the circled 
point at about K=12.4 is an outlier in Fig[2]and may be a false 
measurement due to stellar crowding. Bottom: Stars in the PCI 
region with at least two data points. Nine known variable stars 
(tabulated in Table|2jl are marked with open circles. 



2.2. Deriving <K> 

The epochs of individual K-band frames are distributed over 
20-50% of RR Lyrae star light curve cycles. Uncorrected, 
the measured average magnitude for each variable might dif- 
fer from the intrinsic average value which is required for the 
PLk relation. This fact is also true for the deep/composite 
frames that are considered here for signal-to-noise reasons. The 
adopted way to tackle this issue is taken from Jones et al. (1996; 
referred to hereafter as J96): based on template near-IR light 
curves for both RRab- and RRc-type variables, the method al- 
lows the offset from the true mean magnitude to be estimated 
if the period, optical amplitude, and light curve phase of each 
exposure is known. By taking frames from different light curve 
phases, template shifting acts to minimize deviations from the 
true mean magnitude arising from uncertainties in emphemeris 
phase (J96). For this reason, the magnitude offset associated 
with each composite frames is determined by calculating the 
offset associated with each of its constituent frames, and then 
summing these values, weighted by the constituent integration 



Fig. 4. K-band zero point offset assessment. (V555, V555- 
Kpci) CMD (squares) where Kpci magnitudes have been 
shifted to match the lower end of the RGB fiducial ridgeline 
(dot-dashed/thick) deduced from the work of Ferraro et al. 
(2002; referred to hereafter as F00) taking a distance modulus 
of 15.0. The Cohen et al. (1978; referred to hereafter as C78) 
giant branch data (asterisks and diamonds) has been considered 
for the purpose of zero point validation; taking the C78 data as 
comprising asymptotic giant branch (AGB) stars at V > 13.5, 
the 2nd order polynomial/least squares fit to the C78 data (ex- 
cluding the outliers (diamonds) falls within ± 0.05 mag of the 
apparent AGB stars in the Kpci data. A ±0.1 K-band magni- 
tude envelope (dotted) has been included for convenience. See 
the final paragraph of Sect. l2.1l for more information. 



times. The magnitude offset associated with each star in any 
of the deep frames is in the range -0.05 to 0.03 mag. The er- 
ror reported for the J96 recipe is typically less than 0.03 mag; 
however, for each deep frame there are eleven frames covering 
a range of light curve phases, leading to a negligible average 
error, and is ignored. 

Zero-pointed photometry for a sample of seven RR Lyrae 
stars, labelled in Fig. ^according to Bakos et al. (2000), was 
matched to their associated periods (see Tabled. The RR Lyrae 
stars comprise five stars of the type RRab and two stars of the 
type RRc. Period data with a precision better than 1 % is taken 
from CC01, Corwin (2002, in preparation) and Strader et al. 
(2002). The epochs in Table |2] are taken from CC02 which is 
based on their earlier work (CC01). The CC01 study made use 
of observations made on six nights in 1992, seven nights in 
1993 and one night in 1997. The observations for the present 
study were made in 2000. Of the inner region variables consid- 
ered in this paper all have a precise emphemeris, (except V253 
for which none has been reported). It turns out that the results 
in this paper are insensitive to a random choice of emphemeris 
for this star, as might be expected for an intrinsically low ampli- 
tude variable. No period change has been measured by CC01 
for the inner region sample considered here, although period 
change rates up to about 0.5 days Myr -1 were found for 38 
outer region RR Lyrae stars. Even if such changes were com- 
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mon to each star, it is reasonable to believe that the ephemerides 
that CC01 published are valid for 2000. Incidentally, the ro- 
bustness of the PLk technique to period uncertainties is further 
strengthened by the template shifting mentioned earlier which 
causes considerable reduction in errors arising from errors in 
the light curve phases. B-band amplitudes with a precision bet- 
ter than 5% are available from CC01 for two of the variables 
(A(B)vi93= 1.6; A(B)v20i= 1-6) and has been estimated for 
the remaining five variable stars using the A(B)-period plot of 
CC01 . As noted later in Sect.[3] adopting a conservative error of 
±0.3 mag for any of the B-band amplitudes, produces a change 
in the measured slope of the period-luminosity relation of up 
to ±0.01 mag and an even smaller error in the zero point of the 
relation. 

Another issue is the possible presence of Blazko variables 
in the inner region sample of RR Lyrae stars - the light curve 
shapes of such Blazko-type stars could change on timescales of 
weeks or months. As none of the stars are labelled as Blazhko- 
type variables in CC01 or CC02, there is no 'a priori' reason to 
suspect the presence of such variables in the available sample 
but the possibility remains. However, inspection of CCOl's B- 
band light curves for Blazhko variables indicates that typical 
amplitude changes are up to about 0.5 mag, in which case the 
derived PLk would be relatively insensitive to the effect. 

Next, the intensity weighted average of the magnitudes de- 
termined from the three deep frames is considered. For stars 
with at least two valid measurements, the adopted recipe for 
intensity weighted averaging of photometry is this: the magni- 
tudes for each star are converted into a flux level with a cer- 
tain zero point, averaged and converted back to magnitudes. 
The rms value was calculated for each star; they vary in the 
range 0.02 - 0.05 mag. For line fitting later, the instrumental 
magnitude error for each RR Lyrae star is taken as the average 
(arithmetic mean) of the rms values for the seven variable stars, 
scaled down by \/3, leading to 0.03 mag. In this way, statisti- 
cal differences in magnitude errors are smoothed out. The rms 
value is included in the measurement of the slope of a linear 
fit to the period versus K-band relation described in the next 
section. 

The error in the absolute zero point of the K-band pho- 
tometry from three deep frames is ±y3x 0.01 mag ( relative 
zero pointing error); plus the error in the zero point offset is 
±0.05 mag (systematic), leading to an absolute zero point error 
of ± 0.053 mag. 

2.3. Results 
2.3.1 . Inner region 

The fundamental period P D has been calculated for first- 
overtone (RRc) periods using the conversion LogP Q = 
Log (Prrc) + 0. 127 (Iben 1974), and is plotted against the cal- 
culated mean K-band magnitude for each star in Fig. [5] (top- 
left). Of the nine RR Lyrae stars with K-band magnitudes two 
stars VI 92, and V262 are outside the ± la envelope in Fig.|2| 
as they are systematically too bright; only V262 is plotted in 
Fig.|5]for the sake of presentation but the measured mean mag- 
nitude and uncertainty are given for both in Table [2] VI 92 is 



displaced to a brighter magnitude than the longest period RR 
Lyrae star; CC01 classified the star as an RRab-type variable 
with P = 0.481d, B = 14.94, (B-V) = 0.554, and noted it as a 
blend in their data. For comparison, RR Lyrae stars observed in 
M3 by CC01 have B-band magnitudes as bright as ~ 15.8. As 
V192 is an outlier it is ignored in subsequent discussions. The 
second outlier is V262 which is marked with an open square 
in the left panel and has a small standard error (0.02 mag). For 
the seven remaining RR Lyrae stars there are three mean mag- 
nitude estimates, and the intensity weighted average is taken; 
one exception is V221 for which one measurement appears to 
be spuriously faint by 0.3 mag and the remaining two are used 
to estimate the average. For the measurement of the slope of the 
PLk relation, ctk, is taken to be ±0.03 mag (rms internal scat- 
ter). Both the data and the line fit are shown in Fig.|5](left-top). 
Coefficients of the relation are tabulated in Table[5] The error in 
the true zero point of the relation is made up of the error in the 
absolute zero point (±0.053 mag) of the K-band photometry 
plus (in quadrature) the error in the PLk zero point calculated 
by the line fitting (±0.05 mag) which leads to ±0.073 mag. 
An additional error must be included for the MK-period re- 
lation. For reasons explained in Sect. 13. II the adopted distance 
modulus is 15.07±0.07mag, which leads to a zero-point error 
of <TM K ,zcropt.=±0.10mag in the relation. For different future 
distance modulii this uncertainty must be changed accordingly. 

Out of academic interest, the relation for the photomet- 
rically superior DG data set (FWHM — 0.14") is plotted in 
Fig- E (top-right), shifted by -0.1 mag±0.05mag to match the 
Kpci -band photometry (zero pointed to C78); conservatively, 
the uncertainty in the DG photometry is taken to be 0.05 mag 
for the line fitting. The outlier marked with an open square at 
K=15.87 is V215 in CC01 and has been ignored during line fit- 
ting. There is a significant difference between the best fit line in 
this panel and that of the outer region relation (bottom), most 
likely because no J96 correction has been applied as the ob- 
serving epoch is absent. 

2.3.2. Outer region 

The work by Corwin & Carney (2001; referred to here as 
CC01) provides up to date information on RR Lyrae star 
periods and classification for the outer cluster region. The 
Longmore et al. (1990; hereafter referred to as L90) data was 
reanalyzed using the new periods for the cases of inclusion and 
exclusion of stars classified as Blazhko or double-mode (RRd) 
RR Lyrae stars. In this way, the sensitivity of the slope to these 
stars can be estimated, in addition to improving the precision 
of the relation. K-band magnitudes were taken from L90 which 
have an error ctk = ±0.06 mag, and the latest period values 
come from CC01. There are 48 stars with K-band magnitudes 
from the L90 study, and they cover the radial range 2.6-13.3 arc 
minutes. Of these stars, 37 are now known to be of type RRab, 
nine are of type RRc, one is of type RRd and one, VI 13, is not 
in the CC01 catalogue of variable stars in M3. However, VI 13 
has been recorded by Szeidl (1965) as an RR Lyrae star and is 
therefore included. The period of VI 13 used by L90 is 0.513d; 
the same value is recorded in the Clement (1997) catalogue of 
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variable stars. VI 13 is not classified by Clement (1997) but 
is fairly safe to take it as an RRab-type variable based on its 
period. The RRd star, and 15 of the RRab stars which are clas- 
sified as Blazhko (irregular) RR Lyrae stars by CC01 have been 
ignored for the fitting; they are marked by open circles in Fig.|5] 
(bottom). The gradient and zero point of the PLk relation are 
tabulated in Table [3] for the cases of inclusion and exclusion 
of irregular RR Lyrae stars (rows 2 & 3 respectively). The er- 
ror in the zero point of the absolute PLk relation consists of ± 

0. 06 mag (internal zero point error); ± 0.04 mag (rms fitting er- 
ror); and ± 0.07 mag (distance modulus uncertainty - referred 
to in the previous section). This leads to a zero point error of 
±0.10 mag. 

3. Discussion 

3. 1. Harmony in the PLk relation 

Together with accurate period data from the literature, the K- 
band observations prove on good photometric grounds that 
there is no apparent difference between the inner and outer re- 
gion PLk relation. This conclusion is based on the following: 

1 . The gradient of the inner region relation is weakly depen- 
dent on the way <K> is calculated (intensity- or magnitude- 
weighted average) as shown in Table [3] The agreement is de- 
pendent on the J96 magnitude correction for RRab-type stars; 
excluding the correction produces a significantly steeper slope. 
The agreement between the inner and outer region relation is 
not dependent on uncertainty in A(B); an error of ±30% in any 
of the amplitudes produces a change in slope of ±0.01, which 
is negligible compared to fitting errors, and an even smaller 
error in the zero point. For both the inner and outer region rela- 
tions, periods from CC01/CC02 (and Clement 1997 in the case 
of VI 13) are well determined to a precision significantly better 
than 1%. CC01 mean B-band magnitudes are precise to at least 
5%. 

2. The absolute K-band extinction, Ak, can be derived using 
A K /E(B-V) = 0.13 (Cardelli et al., 1989). As E(B-V) = 0.01 
(Dutra & Bica, 2000) which is negligible, no correction has 
been made for K-band atmospheric extinction which would be 
of the order of 0.001 mag. 

From these findings I conclude the following: 

1. Including the irregular variables, but taking improved pe- 
riods, there is good agreement between the L90 slope mea- 
surement -2.35±0.15 for the outer region and value from an 
updated analysis of the L90 data (-2.35±0.14). There is also 
a fine agreement with the ensemble average of the relation - 
2.338±0.067 determined by Frolov & Samus (1998) who re- 
analyzed the period/K-band magnitude data tabulated by L90 
for nine galactic globular clusters. 

2. Inclusion or exclusion of the irregular variable stars has 
marginal impact on the gradient and zero point of the relation 
(see Table|3] rows 2 & 3); this arises because the irregular stars 
are well-distributed about the PLk line. 

3. There is good agreement between the PLk relation for the 
inner and outer regions (see Table|3] rows 1, 2, & 3). 

4. For the absolute K-band magnitude - period relation, the 
distance modulus to taken to be 15.0±0.07mag (Caputo et al. 
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2000, hereafter referred to as COO) and the reasoning behind the 
value chosen is explained as follows: It is true that the distance 
modulus for M3 is uncertain with values in the range 14.8-15.2 
appearing in the literature (COO). Setting observational errors 
aside for the sake of clarity, the uncertainty arises because in 
order to measure a distance modulus in a certain bandpass, the 
average magnitude of a sample of RR Lyrae stars in that band- 
pass, corrected for reddening, is measured through observa- 
tions while the absolute magnitude is estimated through mod- 
eling. Thus, as a metallicity and a luminosity must be assumed 
for modeling, affecting the model stars' magnitude differently 
in different bandpasses, there may be a mismatch between the 
distance modulus derived for one bandpass and the value de- 
rived for another bandpass. However, Bono et al. (2001) re- 
port a way to find an average luminosity for RR Lyrae stars 
by determining the luminosity that gives the same visual and 
IR distance modulus, thereby yielding a mean distance modu- 
lus. For M3 they find (m-M) K = 15.03±0.07mag which they 
note as being in good agreement with the mean V-band dis- 
tance modulus of 15.00±0.07mag derived by COO. As Bono et 
al. (2001) concluded for this point, this harmony proves the in- 
ternal consistency of the pulsation modeling and reinforces the 
belief in the literature that the reddening towards M3 is almost 
zero. Adopting the COO distance modulus for the inner region 
relation, the absolute K-band magnitude for the inner region is: 

M K = -0.96 ±0.10 -2.42 ±0.16 Log P (2) 

There is good agreement with the slope of the empirical 
metal-free relation available in the literature provided by Jones 
et al. (1992), 

M K = -0.88 ± 0.06 - 2. 33 ±0.20 Log P (3) 

which differs significantly from the slope of 

M K = -1.07 ±0.10- 2.95 ±0.10 Log P (4) 

by Skillen et al. (1993). For the outer region, excluding ir- 
regular variables, one has 

M K = -1.07 ± 0.10 - 2.38 ± 0.15 LogP Q . (5) 
3.2. Comparison with theory 

Bono et al. (2001) have recently calculated the theoretical PLk 
relation using non-linear convective stellar pulsation models 
over the range of stellar temperatures covered by the instability 
region. They considered a range of masses, chemical composi- 
tion, and luminosities, and then assessed the sensitivity of the 
relation to these parameters. They derived that absolute K-band 
magnitude of RR Lyrae stars with metal content in the range 
0.0001 < Z < 0.006 is correlated with period and metallicity 
as 

M K = 0.139 - 2.07(LogP o + 0.30) + 0.167LogZ (6) 

with a total intrinsic dispersion of <tk = 0.037 mag, includ- 
ing a luminosity uncertainty of A Log L ~ ±0.04, and mass 
variations of 4%. 
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The empirical study of L90 found good agreement between 
the PLk relations measured for several galactic globular clus- 
ters which suggests that both the evolutionary (i.e. luminos- 
ity) effects and the spread of stellar masses inside the insta- 
bility strip marginally affect that relation. On the theoretical 
side, Bono et al. (2001) concluded that for a certain metallicity 
the predicted PLk relation is marginally dependent on stellar 
mass uncertainties. A mass difference would appear as a zero 
point offset in the PLk relation. In this context, it is worth re- 
marking that CC01 tentatively suggested a possible scenario to 
explain the presence of some optically sub-luminous RR Lyrae 
stars observed toward the central region of M3, namely that 
collisions due to stellar crowding might prematurely halt core 
helium burning, leading to a lower luminosity than normal. 

Another relevant issue is metallicity. Reviewing the litera- 
ture, one has 

M K = -0.647- 1.72 Log P D + 0.04 [Fc/H] (7) 
M K = -0.76 - 2.257 Log P + 0.08 [Fc/H] (8) 
and 

M K = -0.72±0.11-2.03±0.11LogP o +0.06±0.08 [Fe/H](9) 

by L90, Liu & Janes (1990), and Jones et al. (1992) respec- 
tively. Setting [Fe/H] = -1.34dex (Carretta & Gratton 1997) in 
Eq.(|7J and (EJl gives Mk = -0.18 and -0.19 mag respectively 
for Log P Q = -0.3. Eq. ([9} gives a similar result. In contrast, the 
IR flux method that gives Eq. (13 yields Mk = -0.36, which is 
significanty brighter; this is the same discrepancy noted by L90 
and Bono et al. (2001). Next, adopting a distance modulus of 
15.0 mag, tracks of constant Z were derived from Eq. ([6} and 
overlaid in Fig. [5] (bottom). Measurement errors, especially in 
the zero point forbid a useful fit, but there is fair qualitative 
agreement between the data and log Z = -3.04 dex which corre- 
sponds to the adopted M3 metal abundance [Fe/H] = -1.34 dex. 
It is hard to draw a more substantial conclusion from the graph 
due to L90 measurement uncertainties (±0.06 mag) and the 
distance modulus uncertainty (±0.07 mag). Solving Eq. ( |6j 
yields Mk = -0.367 mag in fine agreement with the IR flux 
method (Eq. (|5})- Hence, the PLk relation, a metal-free rela- 
tion, is in good agreement with the predicted PLZk relation, 
meaning that the observations are fairly well matched by the 
results of stellar evolutionary models. This study cannot ad- 
dress the issue of intrinsic dispersion in the PLk relation for 
signal-to-noise reasons. Incidentally, the signal-to-noise ratio 
of stars in the central regions of globular clusters is affected by 
stellar crowding quite severely. To explore the physics of stellar 
evolution and pulsation theory through the PLk relation in the 
central regions of globular clusters, very high angular resolu- 
tion imaging is required, and is becoming available with adap- 
tive optics assisted near-IR imagers, such as NAOS-CONICA 
on the ESO VLT, and later with NIRCam on the James Webb 
Space Telescope. 
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ID 


Epoch 


Period 


Type 


<K> 


A(K) 


A(B) 




(JD) 


(days) 




(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


vl92 


2449090.225 


0.49790 


? 


12.44 b 


_ 


_ 


vl93 


2449090.327 


0.74784 


RRab 


14.35 


0.34±0.01 


1.6 d ±0.1 


v201 


2449090.133 


0.59408 


RRab 


14.73 


0.34±0.01 


1.6 d ±0.1 


v221 


2449090.125 


0.378752 


RRc 


14.77 


0.27±0.03 


0.75± 0.3 


v241 


2449090.014 


0.59408 


RRab 


14.58 


0.28±0.03 


1.0 ± 0.3 


v252 


2449090.157 


0.50155 


RRab 


14.64 


0.33±0.03 


1.5 ± 0.3 


v253 




0.33161 


RRc 


14.89 


0.27±0.03 


0.75± 0.3 


v262 




0.5647 a 




14.41 c 


0.29±0.03 


1.1 ±0.3 


v264 


2448755.093 


0.35649 


RRc 


14.88 


0.27±0.03 


0.75± 0.3 



Inner region variable stars. 

Cols. (1) ID. V241 and V252 are X17 and KG4 respectively in CC01/CC02; (2) Julian date from CC02; (3) Periods. All from CC02 except 
a Strader et al. (2002); (4) RR Lyrae type from CC02; (5) Intensity-weighted K-band magnitude, ct k = ±0.03 mag (random), plus ±0.05 
(systematic) or ±0.058 mag combined; b ignored as it is an outlier. No J96 correction applied as the magnitude is inconsistent with an RR 
Lyrae star; c is an ignored outlier in Fig.|5|(top-left). (6) Expected K-band amplitudes determined by applying Eq. (7) of J96.; (7) d from CC01. 
A conservative error has been adopted. Other A(B) values were taken with conservative errors from the CC01 A(B) - period graph, their Fig. 9. 
Table 2. Relevant details for nine variable stars from CC02 in the Kpci field. 



Region 


Gradientint 
a 


Zero pointint 
b 


Gradient mag 
a 


Zero point mag 
b 


Comment 


Inner 
Outer 
Outer 
Outer 


-2.42 (0.16) 


14.04 (0.05) 


-2.40 (0.17) 
-2.38 (0.15) 
-2.35 (0.14) 
-2.35 (0.15) 


14.04 (0.05) 
13.93 (0.04) 
13.95 (0.04) 
13.95 (0.04) 


this study 

new analysis of L90 data, excluding Blazhko variables 
new analysis of L90 data, including Blazhko variables 
L90 included Blazhko variables 



Table 3. Coefficients of the general relation <K> = a Log P + b. Values in parentheses are the uncertainties; int and mag indicate 
where intensity- and magnitude-weighted average magnitudes respectively have been used. 
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14.0 



14.5 



15.0 



PCI region 
K = -2.42 (0.16) Log Po + 14.040 




14 - PCI region: DG data 

• K = -2.02 (0.07) Log Po + 14.04 (0.18) 

14.5- 
15.0 ; 




15.5 - 



16.0t 



□ 



0.4 



Fundamental period (d) 



Fundamen;al period (d) 



r > 50 arc sec 

K = -2.38 (0.15) Log Po + 13.93 (0.04) 



Z = 0.0015 




1 5.5 L 



0.4 



Fundamental period (d) 



Fig. 5. Linear least squares fit to the K-band magnitude / period data (solid). The dotted lines form the la envelope of the linear 
least squares fit to the data. Coefficient errors of the general relation aLogP + b, arising from the line fitting are given in 
parentheses. Top-left: Diamonds and triangles mark the RRab- and RRc-type RR Lyrae stars respectively. Top-right: DG data 
shifted by -0.1 mag to match the Kpci data (zero pointed to match C78), with no J96 offset applied as the epoch is unavailable. 
The open square marks an outlying star; it is V215 (K=15.87mag in the DG data and is not in the Kpci data). Bottom: Outer 
region data. Z = 0.0003, 0.0007, 0.0011 and 0.0015 tracks are overlaid (from upper to lower tracks). Data points marked with 
open circles are the Blazhko (irregular) stars and one RRd-type RR Lyrae star. See text in Sect. l2.3l for further information. 



